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SECTION 1 
INTRODUCTION 


Time domain ref lectometry measurements in plasmas are 
of great interest. This is particularly true in studying the 
high density plasma layer which forms around a spacecraft as 
it reenters the Earth's atmosphere. It is well established 
that heat transfer rates for spacecraft in reentry can be 
greatly influenced by flowfield ionization levels. 

Knowledge of the ionization level is, therefore, important 
during reentry as a function of body station, especially in 
the forebody region. Current analytical models are 
unverified and incumbered with unproven simplifying 
assumptions and, as a result, fall 3hort of providing 
dependable electron density predictions. 

The Microwave Ref lactometer Ionization Sensor (MR IS) 
Experiment ^ will make use of the microwave reflection 
properties of ionized gases to determine levels of electron 
density in the forebody shock layer of a reentry vehicle and 
to measure the distance from the vehicle at which these 
electron densities occur. Specifically, the experiment will 
locate the onset and presence of critical electron densities 
corresponding to selected carrier frequencies by measuring 
the amplitude and phase of the reflected signal. The 
critical electron density is the density at which a rapid 
increase in reflection coefficient occurs. The distance to 
the critical density point will be obtained through time 
domain ref lectometry using 64 equally spaced transmitted 
frequencies, occupying a total bandwidth of 2 GHz. At the 
receiver, 64 frequency responses are collected and an inverse 

1 NASA experiment originally proposed for 1996 


fast Fourier transform then provides a time domain impulse 
response for the plasma. The propagation delays of the 
signals reflected from the ionized flow field will be 
measured to determine the stand-off distances to the 
location of the critical electron densities detected. The 
stand-off distance information for the vehicle will be 
compared to the data predicted from the computational fluid 
dynamics models for the spacecraft's trajectory. Such a 
procedure should work well so long as the properties of the 
plasma layer are constant over the time required for the 64 
measurements. Since the positions of the critical electron 
densities will fluctuate during measurement periods, it is 
desirable to study the effects of any such fluctuations on 
the time domain ref lectometry response. The purpose of this 
paper is to study the effects of time dependent-fluctuations 
on time domain ref lectometry for a one-dimensional plasma 
sheath. 

The study of plane waves incident on stable plasma 
layers has been investigated by many authors [1-5]. Their 
work has been focused on solving the differential equations 
which describe the electromagnetic behavior of fields within 
these fixed layers. This paper delves into the 
electromagnetic behavior of fields within a time-variant 
plasma layer. Solutions for continuous inhomogeneous layers 
require, in general, numerical techniques. In this paper 
the plasma is represented by a one- dimensional dielectric 
constant which is allowed to vary in the direction of 
propagation. Practical solutions have been obtained using 
numerical methods for solving the Helmholtz wave equation 
describing propagation through the plasma layer. Sequences 
of density profiles have been generated to emulate the 
time-dependent behavior of moving density fluctuations. The 
plasma has been studied for a range of constituent 
transmitted frequencies (64 stepped frequencies), consistent 
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with the MRIS experiment, and the resulting responses have 
been used to synthesize effective time domain responses. 

Broadly speaking, this paper is divided into four 

parts. Section 2 deals wi th t he models used for a given 

plasma layer which can be represented by a dielectric medium 

with a dielectric constant which can have a negative 

imaginary part representing losses as explained in reference 

6. Section 3 concentrate s on r eflection coefficients for 

different plasma models. To simulate the time domain 

ref lectometry response of a fixed one-dimensional cold 

plasma sheath, solutions for the electric field of a 

normally incident plane wave in a specified electron density 

are used. Illustrative cases with simplified profiles with 

exact solutions are presented for the reader in this 

section. Section 4 looks at simulating the time domain 

response of the plasma by mathematically transforming 

reflection coefficient measurements made in the frequency 

domain. From the responses , propagation delays of the 

signals reflected were measured to determine the distances 

to the location of the critical electron densities. Section 

5 investigates electron density fluctuations and their 

effect on time domain responses. For the sake of clarity, 

many of the problems treated are greatly simplified. For 

instance in Section 3, a homogeneous layer and a linear 

layer are examined as illustrative aids for understanding 

the solutions for a inhomogeneous plasma layer. The appendix 

reviews the time domain theory used in this paper. A 

computer program was developed to solve inhomogeneous plasma 

layer problems. For propagation studies, a time harmonic 

- 

wave traveling in the positive z-direction with an e time 

convention is assumed throughout the paper. Normal 
incidence was chosen to simplify problems and to help focus 
on the intent of the paper. Issues and considerations 
regarding the MR IS distance-measuring scheme are discussed 
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in the conclusion. 



SECTION 2 


ELECTRON PLASMA MODEL 


2.1 Wave Propagation 


In the solution of any electromagnetic problem, 
Maxwell's equations must be satisfied. As shown in 
reference 7, these equations can be used to obtain a 
differential equation describing wave propagation through an 
arbitrary medium. This dif fer ential equation known as the 
one-dimensional Helmholtz wave equation is written as 


^E(z) 

2 


+ LTfJ£( z)E(z) = 0 




( 2 . 1 ) 


where E is the electric field intensity, « is the radian 
frequency, m is the permeability of the medium, and £ is the 
dielectric constant of the medium which is shown as a 
function of z, the axis along which the electromagnetic 
fields propagate. For the case where £ does not depend on z 
(purely homogeneous dielectric), a time-harmonic wave moving 
in the positive z-direction, with an e time convention, 

results in a solution writt en a s 

E(z) = E^e ikz + E 2 e“ ikz , (2.2) 

where E^ and Eg are undetermined constants independent of z 
and k, the wavenumber, is 


k = 





)J£ 


£ 

r o 


(2.3) 


£ 

r 


is the relative dielectric constant of the medium and £ q 
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is the dielectric constant of free space. In the solution 
of the wave equation for E as given by equation (2.2), the 
first term represents a wave with magnitude E^ traveling in 
the positive z direction, and the second term represents a 
wave with magnitude E ^ traveling in the negative z 
direction. 

Now consider a inhomogeneous dielectric, where k varies 
with distance such that 


k = 



£ 


r 


(z) 


(2.4) 


and equation (2.2) is no longer valid. The solution cannot 
be conveniently expressed as a forward and backward 
traveling wave as for the homogeneous dielectric. As shown 
later in Section 3, except for certain special cases such as 
a linear dependence of c with z (see section 3.2), to find a 
solution for the electric field in an inhomogeneous 
dielectric, equation (2.1) must be solved numerically. 

2.2 Properties of Plasma 


Let us consider the particular properties of a 
partially ionized, but electrically neutral, gas insofar as 
they affect the propagation of electromagnetic waves. 
Electrons have a natural frequency of oscillation called 
the plasma electron frequency or more commonly, the plasma 
frequency. It is represented by and is defined as 
follows [8]: 


*2 = -VI 

N ^rn 
o e 


(2.5) 


where c is the permittivity of a vacuum, 
o 


m is the mass of 
e 
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an electron, N is the electron density, and e is the charge 

- 6 _ _ 

of an electron. 

If temperature effects are not important, a plasma can 
be modeled as a dielectric as shown in reference 6, where 



( 2 . 6 ) 


<»> is the transmitted frequency. In the general case can 
be a function of position, thereby, representing a 
inhomogeneous plasma. The refractive index of the ionized 
medium is given by [5] 




(2.7) 


Two cases are possible: 

Case 1 


go > go 1 - — >0 n is real 

N GO 


In this case, wave propagation takes place. 


~ ^ Ca se 2 



go < go m 1 - — <0 n is imaginary 

N GO 

In this case, the fields are evanescent and no wave 
propagation occurs. For a signal to be transmitted through a 
plasma it is, therefore, necessary that the frequency of the 
microwave signal w be higher than the plasma frequency 
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A wave, therefore, may propagate into a medium having 

increasing but as approaches <*>, ^approaches zero and 

a criterion for reflection is met. The point at which the 

electron density level causes the plasma frequency to 

equal <o is termed the turning point in this paper. To 

further illustrate the definition of turning point, in 

figure 1 the number N g of free electrons per unit volume 

increases slowly in magnitude, reaches a maximum, and then 

falls abruptly with further increase in distance. A wave of 

a given frequency w would enter the plasma without 

reflection because of the slow change in N . When the 

0 

density N q is large enough, however, = ° > - Then the 

dielectric constants in equation (2.6) vanish and the wave 
is reflected. In figure 1, h^ is the location of the 
turning point. 


2.3 Reflections 


In the previous section, the properties of plasma were 

studied, and it was stated that a plasma can be modeled as a 

inhomogeneous dielectric. To understand reflections in a 

plasma, a formula for the reflection coefficient must be 

developed. Before we investigate reflections in the 

inhomogeneous dielectric model for the plasma we should 

examine a homogeneous dielectric. Consider a homogeneous 

dielectric with relative permittivity c , a plane uniform 

r 

wave progressing in the z-direction and having its electric 
vector in the y-direction is completely specified by 
equation (2.2) as 

E y (z) = E ie ikz + E 2 e _ikz (2.8) 
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where k is the wavenumber defined by equation (2.3) as 


= /*Vo 


(2.9) 


The reflection coefficient at a location z q can be defined 
as a complex number [9] 


nz o ) = 


„ - ikz 

E 2 e ° 

E ie ikz o 


( 2 . 10 ) 


With equation (2.10), equation (2.8) can be written about 

the point z^ as 
o 


E y (z) = E i ^e ik(z Z o> + r(z o )e~ ik(z Z o ) J , (2.11) 

where E i = Eje ikz o is the incident field for the traveling 

wave at z . Taking the derivative with respect to z, 
o 

equation (2. 11) becomes 


*E (z) 

— = E l [ike lk(z ‘ i! o ) -ikr(z o ) e - ik(z - 2 o ) ] . (2. 12) 


We define a quantity p at the point z Q as 
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(2. 13) 


P = 


<?E (z ) 
y o 

dz 


E (z ) 

y v o 


This factor p is proportional to the admittance of the wave 

given by H /E , where H is the magnetic field in the x 
x y x 

direction. The resultant equation for the reflection 

coefficient at the location z can be written, in terms of 

o 

P, as 


ik-p 

T(z ) = . (2. 14) 

ik+p 

For values of p where the magnitude of r is zero in equation 
(2. 11) , the wave simply propagates in the +z direction with 
magnitude E^. For values of p where the magnitude of r is 
unity, the wave is reflected and travels in the -z direction 
with magnitude E^. 

Having considered a homogeneous dielectric, we can 
investigate the inhomogeneous dielectric model for the 
plasma. As mentioned in section 2.1, equation (2.8) is not 
valid for inhomogeneous media and reference is made to the 
original Helmholtz wave equation. The Helmholtz wave 
equation is written here as 

(z) 2 

£* — + to^p£(z)E(z) = 0 . (2.15) 

As noted earlier for arbitrary variation of £ with z, it is 
not possible to find unique forward and backward waves in 
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equation (2. 15). However, the kinds of variation of c with 
z are such that there is a region of constant permittivity 
(free space) near the transmitting source. It is only in 
such a region of constant e that equation (2. 14) is actually 
evaluated after having found R y (z) everywhere using equation 
(2. 15). Note that this wave equation has many solutions and 
equation (2.8) is a solution for a homogeneous medium only. 
To apply this equation for an arbitrarily varying 
dielectric, a solution must be found numerically for the 
ratio of the field expressions in equation (2. 13) and 
subsequently for the reflection coefficient in equation 
(2.14) (see section 3.3). When the reflection coefficient 
for the modeled plasma is calculated at several frequencies, 
the frequency response for the plasma is known for that 
frequency range and the frequency response data can then be 
transformed to produce the time, domain response for the 
plasma. 
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SECTION 3 


REFLECTION COEFFICIENTS FOR UNIFORM, LINEAR, 
AND ARBITRARY INHOMOGENEOUS LAYERS 


3. 1 Uniform Dielectric Layer 


Consider a plane wave incident on a plane uniform 
dielectric layer as shown in figure 3. The incident 
electric field which is polarized in the y-direction 
propagates in the z-direction through free space (Region I) 
and is normally incident on the layer (Region II) backed by 
free space (Region III). Normal incidence was chosen to 
simplify the problem. The electron density profile for the 
three regions is shown in figure 4. In free space (Regions 
I and III) the electron density is assumed to be zero, and 
for the dielectric layer (Region II), extending from z-z^ to 
z-z 0 , it is assumed to be N . The relative permittivity of 
the layer can be expressed by equations (2.5) and (2.6) 
as 


£ d = 1 


N e 
o 


c m 
o e 


(3. 1) 


N q is chosen so that > w and therefore < 0. The 
relative permittivity for this geometry is shown in figure 
5. The field in Region III can be written with unity 
magnitude as 

E (z ) = e ik o z ( z > z 2 ) (3.2) 

^ Electron density profile consistent with experimental 
predictions (see figure 2). 
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and 


*E(z) 

dz 


ik e 
o 


ik z 
o 


( z > z 2 ) 


(3.3) 


At z=z g . the field can be written as 

„ / x ik z 0 
E(z 2 ) = e o 2 


(3.4) 


and 


*E(z 2 ) 

dz 


ik e 
o 


lk o z 2 


(3.5) 


In Region II, 



(3.6) 


and the index of refraction, n can be written as [9] 



(3.7) 


where e is the relative permittivity. In general let 


/ 


£ 

r 




(3.8) 
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where and n^ are the real and imaginary part of the index 
of refraction, respectively. Since n is purely imaginary 
for the relative permittivity here, 


k = ik n. 

o 1 


(3.9) 


The field in Region II can be written as 


E(z) = 


C^e 


ikz 


C 2 e 


- ikz 


( z z < z 2 ) 


(3. 10) 


Using equation (3.9), equation (3.10) can be written as 


n . , r% k n . 2 p kn.z . . 

E(z) = Cj^e o i + C 2 e o l z^ z < z^ 


(3. 11) 


and 


*E(z) 

Oz 


, n -kn.z 
-k n.C.e o l 
oil 


, n kn.z 
k n.C 0 e o l 
o i z 


(3. 12) 


Using the boundary conditions, the tangential electric and 
magnetic fields are continuous at the interface z=Z 2 , we can 
equate the field expressions in equations (3.2) and (3.11), 
and equations (3.3) and (3.12), respectively, to solve for 
and Cg (see reference 2). Note that 
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1 - 


1 

in. 


-2n,k z 0 
e i o 2 


1 + 


in , 


will be very small if n^Zg is several free space 
wavelengths. Choosing Zg such that l^n^Zg-z^) » 1 renders 
l C 2 l « |C^|e ^ n i^o z l. Similiarly, in Region I we can equate 
field expressions at z=z^ and the reflection coefficient at 
z=0 can be written as (see reference 2) 


P(z=0) 



(3. 13) 


where 


n , 


Co = 


3 ' 2i 


and 


• ■ - ®i(»is r K k » ni,i,1 'i=r > 


- "i [ r f 1 + 1 1 k n.z.d-rj-) 

- [~ c 2l 1+ nrJ e ° 1 1 in i 


(3. 14) 


♦ c i( i -iH I K k o n ^i (U ^r ) ] • « 3 - 


15) 


Since |Cg | « | C ^ | e ^ n i^o z l, 
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r( z=0) 


1- in 


i 2ik z 1 
— e o 1 


(3. 16) 


1+in^ 


where 


n 


i " 



(3. 17) 


3.2 Linear Layer 


In the previous section we calculated the complex 

reflection coefficient for a uniform dielectric layer. Now 

consider a linear layer as shown in figure 6. A plane wave, 

polarized in the y-direction and traveling in the z-direction 

through free space (Region I, N =0), is incident on the 

layer at z=z^. As the wave progresses through the layer 

(Region II), it encounters greater electron densities. At 

the outer edge of the layer, z = Zg, the electron density N q 

is at its maximum, N q . As the wave leaves the layer, it 

returns to free space (Region III, N =0). The electron 

© 

density profile can be expressed as 


N N z 1 

N = 2L_ z ° 1 

e z 2 - Zl 


z 2 -z 1 


and in air as 


N = 0 
e 


( z 1 < z ^Zg ) 


(3. 18) 


( z^Zj and z^Zg )• (3.19) 
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When the electron density N is large enough, the relative 

© 

dielectric constant £ in equation (3. 1) vanishes and the 

r 

wave is reflected. Let us call this critical electron 
density value N . For the dielectric constant to vanish 


N = - 


c m 
o e 


(3.20) 


By using equation (3. 18) we find this to happen at the 
turning point 


z = z^ + 


u £ ra Az 
o e 

u 2 

N e 
o 


where 


(3.21) 


Az = Zg ~ z ^ (3.22) 

The field in Region III can be written with unity magnitude 
as 


„ , . ik z 
E(z) = e o 


( z*z 2 ) 


(3.23) 


and 


JB(z) 

9z 


ik e 
o 


ik z 
o 


( z^Zg ) • (3.24) 


In Region II the Helmholtz wave equation derived in 
reference 9 becomes the Airy differential equation as shown 
in reference 2. The field expressions can then be written in 
terms of Airy functions (see reference 2) as 
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E (u) = C 1 Ai(u) + C 2 Bi(u) 


(3.25) 


and 


Ci^jAi(u) * C^Bifu) 


where the variable u(z) is defined as 


u(z) = - 


k 3 Az 
o 


n r 


K, 


1 - 


■§KH) 


and 


K 1 = 


N e 
o 


2 _ 

c £ in 
o e 


(3.26) 


(3.27) 


(3.28) 


Using the boundary conditions at z=z^ and z=z 2 > we can 
equate field expressions and solve for the reflection 
coefficient at z=0 (see reference 2). At z=0 we are in free 
space, a region of constant permittivity, the reflection 
coefficient at z=0 can be written as 


( C 1 L 1 + C 2 L 2 > 21k ° Zl 

r( z =0) = — — (3.29) 

C 1 L 3 + C 2 L 4 


where 



(3.30) 


L x = Ai(u(z 1 )) 


IjjAituCz 


1 


)) 



dz 


Lg = Bi(u(Zj)) 


d 


feBi(u( Zl )) 


ik — 
1K o du 


(3.31) 


Lg = AKufz^) + 


SlAKuUj)) 


Ik 


<^Z 

o 


(3.32) 


L 4 = Bi(u(z 1 )) + 


3^i ( u ( Zj)) 

ik 57 
o #u 


(3.33) 


To illustrate this exact solution for the reflection 
coefficient, a particular linear profile was chosen and is 
shown in figure 7. The electron density begins at the 
front interface, z=0, at a value of zero and rises to a 
value of 1 x1(j u electrons per cubic meter at z=14 
centimeters (the exit point). The real and imaginary parts 
of the reflection coefficient at the front interface are 
shown in the right column in table I for 74 to 75 GHz. 


3.3 Plasma Layer 

In the two previous sections the complex reflection 
coefficients for a constant dielectric slab and a linear 
layer were derived exactly. The solution for an arbitrary 
inhomogeneous plasma layer requires, in general, 
numerical techniques. The plasma is represented by a 


19 



scalar, isotropic, and inhomogeneous dielectric constant. 
Consider a plasma layer with a electron density profile as 
shown in figure 2. An incident electric field, which is 
polarized in the y-direction, propagates in free space 
(Region I), in the z-direction, and is normally incident on 
the plasma layer (Region II) and backed by free space (Region 
III) as shown in figure 8. Again, normal incidence was 
chosen to simplify the problem of studying the effects of 
electron density fluctuations on the time-domain 
ref lectometer response for a one-dimensional plasma sheath. 

The Helmholtz wave equation for the electric field in 
a inhomogeneous plasma layer is written as 


(z) ? 

— + k (z ) E (z ) = 0 (3.34) 


where k, the wavenumber, is 


k(z ) 


<Jk> 


/ 


P«(z) 


(3.35) 


The wave propagates through the plasma (Region II) and is 
transmitted to free space (Region III). Since the relative 
dielectric constant in free space equals one, the solution 
of equation (3.34) in Region III is readily found. 

To begin the solution, the electric field in Region III 
may again be written with unity magnitude as 

E y (z) = e ik o z (3.36) 

The derivative needed to define p in equation (2.13) is 
again given by 
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(3.37) 


<»E ( 2 ) 

y 

dz 


ik e 
o 


ik z 
o 


The field at z = d is assumed to be unity and, therefore, can be 
written as 


E (z=d) = e ik o d 

y 


(3.38) 


and 


<9E (z=d) 

y 

dz 


ik e 
o 


ik d 
o 


(3.39) 


Equations (3.38) and (3.39) serve as boundary conditions 
for solving equation (3.34). Dropping the polarization 
notation and making the substitution 


*E(z) 

S(z) = _ (3.40) 

&z 


equation (3.34) becomes 

+ k 2 ( z ) E ( z ) = 0 (3.41) 

Using a fourth-order Runge-Kutta method [10], equations 
(3.40) and (3.41) can be integrated to find the electric 
field and its derivative. Once these are found at z=0, p at 


21 


z=0 becomes 


<9E ( z - 0 ) 

dz 

P = (3.42) 

E (z = 0 ) 


as defined by equation (2. 13) and the reflection coefficient at 
z = 0, in terms of p, is 


ik-p 

r(z = 0) = (3.43) 

ik+p 


One can note that only the ratio of S(z) to E(z) is actually 
needed to compute the reflection coefficient. Thus, the 
magnitude of E(z) can be kept near unity by normalizing the 
solution at the end of each Runge-Kutta step. 


3.4 Verification of Runge-Kutta Method 

In section 3.2, a linear profile was chosen to 
demonstrate an exact method ( Airy-equation) of determining 
the reflection coefficient. A comparison of the numerical 
Runge-Kutta solution and this exact solution is made. The 
particular linear profile chosen is shown in figure 7. The 
real and imaginary parts of the reflection coefficient at 
the front interface for the Runge-Kutta method and the exact 
method are shown in table I for 74 to 75 GHz. The 
Runge-Kutta method compares favorably with the exact 
solution as can be readily seen in table I. In all cases, 
the number of Runge-Kutta steps used was tested so that 
significant changes in computed values were not observed for 
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larger numbers of steps. 


SECTION 4 


TIME DOMAIN MEASUREMENTS 


4.1 Uniform Dielectric Layer 

As stated, the complex reflection coefficient, for a 
range of constituent transmitted microwave frequencies, can 
be used to synthesize the effective time domain response of 
a plasma layer (see Appendix for a review of time domain 
theory). Sixty four frequencies, for a bandwidth of 2 GHz, 
were chosen to be consistent with the MRIS experiment. 
Reflection coefficient data in the frequency domain are 
shifted and transformed to give a baseband time domain 
response using a decimation- in- time fast inverse Fourier 
transform [11]. The resulting time domain response emulates 
a baseband continuous wave signal, where the peak magnitude 
marks the location of the turning point. The frequency 
domain data are windowed by a Kaiser-Bessel window to reduce 
unwanted interference, and it should be noted that the peak 
response of the data is normalized to the reponse of the 
window. To illustrate the technique of using the time 
domain response to locate the turning point, the results of 
section 3. 1 are used for a uniform dielectric layer. 

Consider a plane wave incident on a uniform dielectric 
slab, with -0.5, located 10 cm from the source, as shown 
in figure 5. A negative permittivity was chosen to simulate 
a purely reactive dielectric reflector. The wave propagates 
through free space (Region I) at the speed of light c and is 
normally incident on the slab at z^lO centimeters. The 
slab is 14-cm thick (Region II) and is backed by free space 
(Region III). As stated by equation (2.10), the ratio of 
the amplitude of the field reflected by the dielectric slab 
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to the amplitude of the field incident at z=z ^ is called the 
reflection coefficient r of the uniform dielectric slab. We 
should note by equation (3. 16) that r is the reflection 
coefficient pertaining to z = 0 and takes into account the 
effect of the propagation path from 0 to z^. The effects of 
the back of the slab are attenuated and not seen. The wave 
is incident and reflected at the z-z^ interface, which is 
the turning point for the dielectric layer. The magnitude 
and phase of the reflection coefficient for the uniform 
dielectric layer are shown in figures 9 and 10, 
respectively, for 74 to 76 GHz. The phase plot is relative 
to the phase at 74 GHz. T hes e reflection coefficient data 
were transformed and windowed with the resulting time domain 
response of the layer as shown in figure 11. The response is 
shifted to the left of time t=0 and the shift is the 
propagation delay of the wave. Note that this delay is 
measured on the negative side of the time axis and has 
negative values. Due to the truncated frequency domain, 
ringing is associated with the time domain response. A 
Kaiser- Bessel window [10] was chosen and the frequency data 
were multiplied by the real Kaiser-Bessel weights to emulate 
a pulse train with very low sidelobes or "ringing" between 
pulses. Sidelobes can limit the dynamic range of the time 
domain measurement by hiding low-level responses within the 
sidelobes of the higher level responses. 

To estimate the location of the turning point, the time 
for the wave to travel to the turning point had to be 
determined. The time domain response was used to calculate 
this propagation time. In figure 12, we see the transformed 
reflection coefficient for -2 to 2 nanoseconds plotted with 
the effective transmitted pulse (i.e. , the transformed 
Kaiser-Bessel window) as a reference centered at t=0. The 
time shift between the two plots corresponds to the 
round-trip propagation delay for the wave and is labeled 
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2t,, where t, is the one-way propagation delay measured on 
d d 

the negative side of the time axis. Since the wave 
traveled at the speed of light c in Region I and t^ was 
measured to be .3352 nanoseconds, the wave traveled 10 cm 
before it was reflected. This is in agreement with the 
geometry of the problem, where the turning point is located 
at z-z i (z 10 cm) . 

4.2 Plasma Layer 

For the plasma layer the turning point is located 
within the layer. To estimate the turning point distance, 
the velocity through the plasma must be obtained. This 
suggests a problem since the plasma is modeled as a 
inhomogeneous dielectric and the velocity of the propagating 
wave is dependent on the permittivity of the medium. Some 
means must be found to approximate the velocity profile. A 
possible solution is to assume that regions of the plasma 
profile are piecewise linear to find average velocities for 
each such region. To approach this problem, let us look at 
the linear layer discussed in section 3.2. The electron 
density profile of the layer bounded on both sides by free 
space is shown in figure 6. The phase velocity through the 
layer is 


/ \ c 

v p (z) = — 

T^ r (z) 


(4.35) 


The group velocity which describes the transport of energy 
is related to the phase velocity and is 
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At the front interface of the layer (z=z 1 > the relative 
permittivity is that of free space and is equal to one. The 
group velocity there is simply the speed of light c. As the 
wave travels through the layer , it encounters greater 
electron density N values and its group velocity decreases. 
When the electron density is large enough, the relative 
dielectric constant e vanishes in equation (2.6) and the 
phase velocity becomes infinite forcing the group velocity 
in equation (4. 36) to become zero. Note that this occurs at 
the turning point. It is, therefore, possible to assume that 
the average group velocity, as the wave propagates from the 
front interface to the turning point, is c/2. For 
investigations of wave propagation in a linear layer plasma, 
such an assumption is, therefore, quite permissible. However, 
as the plasma becomes nonlinear, the c/2 assumption becomes 
less accurate. 

To illustrate this point, shown in figure 2 is an 
electron density profile which is consistent with the MFIS 
experiment that is being simulated. As an example, let us 
assume the profile to be linear for electron densities less 
than 1. 00xl0^ 4 /cm 3 . We will assume this to be true only for 
this region in order to derive an average group velocity for 
the region. Note that for electron densities less than 
1. 00xl0 14 /cm 3 the critical frequency is less than 90 GHz by 
equation (2.5). For this region, as before, the average 
group velocity is approximated as c/2. For 74 to 76 GHz, the 
time domain response was obtained. With this result, the 
turning point was estimated at . 77 cm. This is in 
reasonable agreement with the exact location of . 87 cm . It 
should be emphasized that this approach was given here as an 


example and that a better approximation could be made by 
further segmenting the same region to find more accurate 
velocity values. Such further segmenting would require 
additional measurements at frequencies below 74 GHz. 
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SECTION 5 


PLASMA FLUCTUATIONS 


5. 1 Background 

In the previous chapter we synthesized the effective 
time domain ref lectometry response of a plasma using a range 
of transmitted microwave frequencies. These results as a 
whole gave us an estimate for the location of the critical 
electron densities in static plasma models. However, 
macroscopic effects of electron density fluctuations, which 
are not included in standard aerothermodynamic simulations, 
may have a noticeable effect on hypersonic reentry flow 
fields [12]. In order to monitor these effects, we will 
extend the foregoing discussion on the propagation of waves 
through static plasma models. Thus, we consider the 
important effects of time-dependent electron density 
fluctuations on the time domain response of a plasma. 

Fluctuation of the profile results in motion of the 
turning point and is emulated between frequency step 
measurements. Motion occurs between frequency step 
measurements , assuming that each individual measurement is 
accurate at each step. Two Doppler effects resulting in 
motion are studied; the first type of motion involves 
reducing the electron density profile while preserving the 
shape of the profile and the second type of motion involves 
modulating the density profile. The first effect is 
demonstrated in figure 13 where a sample electron density 
profile is reduced, moving the turning point from h^ to 
and then h^. The second effect can be seen in figure 14 
after the same electron density profile is modulated. Two 
different rippled profiles are shown moving the turning 
point from d^^ to dg. Both effects are applied 
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incrementally between frequency steps and consequently 
results in erratic motion of the turning point. The 
nonuniform motion of the turning point is estimated by 
calculating the location of the turning point at specific 
tiroes during the measurement period. These discrete 
calculations give an indication of the overall motion of the 
turning point caused by the two effects. Turning point 
velocities achieved are significantly less than predicted 
flow field Mach velocities of hypersonic reentry vehicles 
[13]. 

5.2 Reduction of Electron Density Profile 

The first type of motion resulted by reducing the given 
electron density profile values of figure 2 by percentages 
without changing the shape of the profile. Shown in figure 
15 are the time domain responses for three electron density 
profiles for 74 to 76 GHz. Time t=0 corresponds to the 
location of the time domain response of the window used. 

The first profile has 100 percent density values, the second 
has 95 percent values and the third has 90 percent values. 

It is evident that the responses shift toward the left, 
moving the turning point left of time t=0, with decreasing 
percentages. To understand the displacement of the turning 
point, the location of the turning point was calculated for 
the 100 percent case and the 90 percent case at 75 GHz. By 
using equations (2.5) and (2.6), the turning point moved 
. 13 cm for the 10 percent variation in the profile. Using 
the dwell time of 2.5 milliseconds for the proposed MRIS 

o 

instrument , a turning point velocity of . 52 meters per 
second was computed. 

O 

MRIS Experiment Requirements Document (ERD) 
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To emulate a moving profile response, where the turning 
point moves as the frequency is stepped from 74 to 76 GHz, 
the profile was reduced from 100 percent to 90 percent. The 
time domain response for this emulated variation in the 
profile is shown in figure 16, left of a. static 95 percent 
profile response. The 95 percent profile response is 
considered to be the response at which the turning point is 
at its average position for the 10 percent variation in the 
profile. A shift between this average and emulated response 
of . 17 nanoseconds was computed. By using an average 
velocity of c/2, the turning point was found to have shifted 
1.28 cm in comparison to the variation in the static profile 
turning point of . 13 cm. This significant increase in 
profile shift due to the media motion is a major potential 
error. 


5.3 Modulation of Electron Density Profile 

For the second type of motion investigated, modulation 
was introduced into the profile. An expression is written 
in the form 

1 + A cos(Kz + v>) (5. 1) 


where A is the amplitude, w is the phase of the modulating 
wave, and K is the spatial wavenumber defined as 



where x is the spatial wavelength. This function was used 
to modulate the profile. By varying the phase w of this 
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function, the disturbance described by equation (5. 1) was 
3et into motion and traveled across the profile. 

A modulating wave with an amplitude of . 05 and a 
wavelength of 2 cm was chosen. To mimic a modulation 
velocity of 2 meters per second, the phase v' of equation 

o o 

(5. 1) was varied from 0 to 90 as the frequency was stepped 
from 74 to 76 GHz. The time domain response is shown in 
figure 17. By squaring the transform reflection coefficient 
in figure 17, the power distribution in the time domain can 
be viewed. The power distribution for the response in 
figure 17 is shown in figure 18. Greater modulation 
velocities up to 8 meters per second were investigated and 
results for the reflection coefficient and the power 
distribution are shown in figures 19 and 20, respectively. 
Further work was done for 140 to 142 GHz with the same 
modulating wave to illustrate greater effects at higher 
frequencies, the transform reflection coefficients were 
determined, and the associated power distributions are 
shown in figures 21 and 22. In addition, for 74 to 76 GHz 
the amplitude of the modulation was raised to .25 and 
results are shown in figures 23 through 25. 

Smearing and shifting of the time domain reconstruction 
is quite evident for the cases studied. For example, at 140 
to 142 GHz the response maximum moved from -1 nanoseconds to 
-2.4 nanoseconds as the modulation velocity was increased 
from .22 meters per second to 3.33 meters per second (see 
figures 21 and 22), and as the velocity was increased further 
to 8 meters per second, the response became more smeared and 
shifted with multiple echoes (see figure 22d). Various 
methods for finding the turning point from these responses 
can be employed. Two approaches are presently realizable; 
one approach involves finding the peak power level and 
designating it as the location of the turning point and 
another approach involves finding the centroid of the 
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area under the power distribution plot, where the centroid 
would mark the location of the turning point. There are 
a few problems with these approaches. Smearing of the time 
domain response results in several peaks which may be used 
to define the turning point, limiting the accuracy of the 
first approach. Shifting of the response centroid would 
also limit the success of the centroid approach. Using the 
previous example for 140 to 142 GHz, the centroid of the 
distribution moves from - 1 nanosecond to -2.6 nanoseconds as 
the modulation velocity is increased from .22 to 3.33 meters 
per second. 

From the results presented, a limitation of 2.67 meters 
per second on the turning point velocity is estimated for 
the accuracy of the instrument at 140 to 142 GHz with 

modulation amplitude A= . 05 . For 74 to 76 GHz with 

modulation amplitude A=.25, it is estimated that there 
is potential error with turning point velocities greater 
than 4.67 meters per second. 


The following conclusions should be kept in mind when 
assessing the location of the turning point: 

1. Amplitude fluctuations in the electron density profile 
cause a shift or delay in the effective time domain 
response. 

2. Smearing of the time domain response, creating multiple 
peaks, becomes significant for rapid fluctuations in the 
electron density profile. 

3. Strategies for locating the turning point can be chosen 
to possibly minimize adverse effects resulting in shifting 
and smearing of the time domain response. 
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SECTION 6 


CONCLUDING REMARKS 


The results of studying the effects of moving electron 
density fluctuations on frequency generated time-domain 
ref lectometry in a one-dimensional plasma layer have been 
presented. On the basis of the studies made and the data 
obtained, the following remarks and conclusions may be made: 

(a) Different models for a plasma layer have been developed 
and discussed in length. 

(b) Equations describing wave propagation through different 
plasma models have been formulated. 

(c) The computer program synthesizes, using 64 frequencies, 
the time domain response of a given plasma electron 
density profile. A linear profile with an exact 
solution (Airy solution) compared accurately with 

the program's numerical Runge-Kutta solution. A 
uniform dielectric slab with a known turning point was 
also used successfully to verify the code. 

(d) The average velocity of an electromagnetic wave 
propagating through plasma must be accurately estimated 
to determine the location of critical electron 
densities. 

(e) Varying electron density levels corrupt time domain 
and distance measurements. In this work it has been 
shown that lowering or reducing the electron density 
levels of a given electron density profile, while 
maintaining the shape of the profile as in figure 13, 
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results in motion of the turning point, and the 
effective motion has a significant effect on measuring 
critical electron density locations. 

(f) Modulating an electron density profile with a waveform 
creates a disturbance or ripple adversely affecting the 
time domain response of a plasma. Waveforms with phase 
variations emulating motion across the profile were 
used and greatly influenced simulated measurements, 
especially for rapid phase variations. 

(g) A technique such as the centroid method for locating 
the turning point may be used to reduce the effects of 
electron density fluctuations on turning point 
estimates . 

All of these issues, and perhaps more which now remain 
unidentified, must be addressed and quantified in order to 
arrive at an estimate of the usefulness of time domain 
ref lectoraetry for locating critical electron densities in 
plasma. 

A 10 percent reduction of the electron density profile 
as described in section 5.2 shifts the turning point 
significantly and may contribute to error. By modulating 
the electron density profile as in section 5.3, there is a 
potential for error when turning point velocities greater 
than 2.67 meters per second are achieved at 140 to 142 GHz 
with modulation amplitude A=. 05. For 74 to 76 GHz with 
modulation amplitude A=.25, turning point velocities greater 
than 4.67 meters per second may also induce errors. It must 
be concluded that a distance- measuring scheme using time 
domain ref lectometry (i.e., MRIS) could become inaccurate if 
some of the plasma fluctuations investigated in this paper 
are encountered. It has been the intent of this paper to 
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help identify and possibly solve similar induced errors 
should they occur. 


Appendix 


Time Domain Theory 


Time domain theory plays an important role in the MRIS 
distance-measuring technique. The response of plasma to a 
time domain signal is used to estimate distances to critical 
electron densities. Plasma reflection information in the 
frequency domain is taken and transformed to the time domain 
where propagation delays are measured to calculate these 
distances. In this process, a frequency down conversion 
takes place so that the resulting time domain response 
emulates a baseband continuous wave signal. Further, an I,Q 
(In-phase, Quadrature) detector is used to construct the 
time domain signal. Processing of this signal is done in 
the form of windowing to reduce unwanted interference. 

To understand exactly how an I,Q detector can obtain a 
time domain response by only using discrete samples of 
amplitude and phase in the frequency domain, consider first 
the problem in reverse. Let a transmitted periodic pulse 
train in the time domain be represented by the signal 


M+k 

f(t) = E a e 
n=M n 


-inw 

o 


t 


-(M+k) 
+ E a 
n=-M 


-in<*> t 
e o 


n 


(A. 1) 


with 

<o q = 2rr/T o (A. 2) 

where is complex and T is the period. The signal is 
real and periodic with 2k spectral components, where k is 
the number of spectral lines in the positive or negative 
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frequency domains. The amplitude spectrum starts with a., in 

M 

the positive frequency domain and a in the negative 
frequency domain. M is chosen here a3 an arbitrary 

9 

constant. Let n=-n in the second summation and drop the 
primes. The signal becomes one summation 


M+k 


f (t) 


= E fa 

n=M n 


- in<*> t 
e o 


mw 

a e o 
-n 




(A. 3) 


Let a = 
-n 


% 

a so that f(t) is real, then 
n 


f(t) = 


M+k r 

= E l 

n=M'- 


% 

(a +a )cos n<*> t 
n n o 


- i (a -a )sin 
n n 


nut 
o , 


(A. 4) 


Further, let 


A n = 


a + a 
n n 


and 


B = i fa - a*l 
n \ n nj 


(A. 5) 


(A. 6) 


To make a frequency down coversion of the signal we 
start by introducing the sinusoidal tone cos poj^t. For any 

p where 

M < p < M+k and k < M 
multiply both sides by cos p^t, thus 


38 


M+k 

f ( t ) cos pw Q t = E [A n (cos nw Q t*cos pw Q t) 
n- M ^ 


B n (sin nw^fcos p^ Q t)J 


M+k , 

E K |cos(n-p)<*> o t + 
n=M 2 


cos (n+p)<*> 




B 


tj- ^sin(n+p)w Q t + sin(n- 


n-p)<o o t]J . 


(A. 7) 


Low pass filtering this response for frequencies less than 
2Mt»> .we obtain a filtered version of f(t)cos pw t as 

O J u 


r ^ M+k 

f^n cos(n-p)<*> t - ®n sin(n-p)<*> tl 

If (t) COS poj 1 1 = E 

l* ° J n=M 

L~2 ° 2 J 


(A. 8) 


where LPF denotes the low-pass filtering. The frequencies 
present range from 0 to kk> Q . 

Now let f(t)cos P w 0 t be placed in a narrow band low 
pass filter such that frequencies of <*> Q and above are 
completely cut off. Then all of the terms in the summation 
except for n=p are suppressed and 

. A 

BPF ff(t)cos P^tJ = ~2 • (A. 9) 

where BPF denotes the bandpass filter. This signal would be 
detected in the I (In-phase) channel of an I,Q synchronous 
detector with the heterodyne frequency being P w 0 » or 
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equivalently in a homodyne synchronous detector if the 
signal f(t) represented a single frequency signal of 
frequency po> Q rather than a pulse train signal. 

To find the output of the Q (Quadrature) channel 
multiply f (t) by sin po^t. Similarly we find 

M+k f 

f(t)sin P w Q t = E Ia^cos nw^fsin p<*> Q t 

n=M 

- B sin nw t*sin p<*> t| 

n o o J 

M+k „ 

n=M |' _ §[ sin ^ n+ P^ < * > 0 t + sin(p-n)« o tj 

- ^n Icos (p-n)<*> t - cos(n+p)<*> t 

2 L o o 

(A. 10) 



After filtering out frequencies less than 2Ma> Q 


LPF |f(t)sin pw^tj = 


M+k 

E 

n=M l 


n sin(p-n)w t 
— o 


B 

~2 cos (p-n)<*> o t 


(A. 11) 


and after cutting off at and above 


BPF |f(t)sin poj Q tj = 


B 

_ _E 
2 


(A. 12) 


The quantity -B p /2 is then the output of the Q channel. If 
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the medium (plasma) being probed is not dependent on time, 

then the sequence of, essentially d.c. measurements, and 

B would be completely sufficient for the construction of 
P 

f (t). 

With the construction of f(t) by the I,Q detector, 
further processing is done to eliminate interference. To 
show how this is done, consider a real function, later to be 
called a "window" function 

M+k 

w(t) = E W cos n<*> t (A. 13) 

n=M n 


and form another function 

T /2 
o 

g(t) = J f(r)w(t-T) dr 

-T /2 
o 

where 

T = for all < t < — 

o co & & 


Let g(t) be periodic with period T q so that 


V 2 r M+k 

g(t) = f ^ E A n cos nw Q T - B n sin nu Q T 


r M+k l 

I r W cos m<*> (t-r) I dr 
«■ m=M ra ° J 


(A. 14) 
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M+k M+k T o /2 

= E E J 

n~ M n5— M _ rp y>2 


I |A W cos nco r • cos m<*> (t-r) 
[inn o o 


( B W sin nw t*cos ra<*> (t-r) dr . 
n n o o JJ 


(A. 15) 


The first terra can be expanded where 


cos nw t • cos raw (t-r) = 
o o 


cos raw 


+ s 


t£ gcos(n-m)w o r + ^cos (n+ra) w^t J 

in raw 1 1 isin(ra+n)w r + isin(ro-n)w r 1 
o l ^ o c o j 


(A. 16) 


Note that integration over a complete period r causes all 
terms to vanish except those where ra=n. Expansion of the 
second terra in equation (A. 15) gives a similar result. 
Thus , 


T M+k 

g(t) = -4 E W 

^ n=M 


multiplication terra- by- terra of the Fourier series 

representing f(t) by W n> produces the Fourier series of g(t) 

given by equation (A. 14). Thus, any pulse train described by 

a set of A 's and B 's and, hence, f(t) can be transformed 
n n 

into another pulse train with a differently shaped pulse by 

multiplying each terra in the Fourier series by W . 

As an example of how a signal f(t) is windowed, let a 

set of A 's and B 's be chosen arbitrarily and let 
n n 

A =l/k for all M^n^M+k and B =0 for all such n. Thus, using 
n n 


■I A cos 
n ^ n 


nw t 
o 


- B n ain 


nw t 
o 


} 


(A. 17) 
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equations (A. 5) and (A. 6) in equation (A. 4) gives 


M+k 

f(t) = E 

n=M 


Note that f(t) is clearly periodic with period 2n/u> Q and has 
equally weighted spectral components. Consider the function 
f(t) for -n/a> o < In order to illustrate the windowed 
function f(t), consider a limiting form for this function 
obtained by letting co q -»0, M-*co, k-»oo in a certain way so that 
the number of components in the spectrum becomes infinite. 

We first write 


M+k ^ M+k « T 

f (t) = E Id cos no) t = E ^Tk 003 n V ’ (A - 19 > 

n=M v J n=M o o 


T M+k 

2 §k Z (cos n V )“ 0 ■ T 0 = I" • (A- 2°) 

n- n o 


so 

f(t) = 


j^lcos 


nto t 
o 


(A. 18) 


Then let <*> -*0, M->oo, k-*» such that 
o 


noo = oo 
o 

(A. 21) 

oo :Au 

(A. 22) 

O 


Mco rO, 
O 1 

(A. 23) 

(M+k)oo o =Q 2 

(A. 24) 

T 

o _ m _ 2 n 

~u - 1 - 
K 0^-0, 

(A. 25) 

2 1 


can write f^(t) for this limiting case 
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v" ■ h r 

n. 


cos <*>t dw 


T 1 
2 


- £ [sin 0 2 t - sin O lt ] 

(A. 26) 


Now 


sin Ci^t ~ 3 


in t = 2|sin| °2 °1 |t j cos — — - 


(A. 27) 


so the limiting form of f(t) is 


f 2 (t) 


_o 

rr 


r ' 

r 

o n -a. 


2 1 


2 J 

k. 


V°1 

sin — k t 


°2-°l 


O +Q 

cos — 7 ) t , (A. 28) 


a carrier at (O^+Ci^) /2 modulated by a sine function whose 
first zero is at t= 2n/(ft 2 -0^). The limiting form of the 
pulse in equation (A. 28) may be repeated such that 
f(t+nT ) z* f(t) for all integer n. In reality this limiting 
form of the pulse shape is not achieved. We must really 
deal with the finite sum and finite of equation (A. 18). 
Thus, the "sine" modulation is only suggestive of what the 
actual pulse shape using finite will be. 

The use of no window implies that we are using a pulse 
train exactly given by equation (A. 18), which in the 
limiting case looks like a "sinc M -modulated carrier pulse 
train given by equation (A. 28). Other shaped pulse trains 
can be produced by windowing f(t), to give 


44 


t 

o 


(A. 29) 


g(t) 


M+k 

E 

n=M 



rcos 

k 


The Kaiser-Bessel set of frequency weights W n gives a well 
defined pulse shape that is nearly zero over roost of the T q 
period except for the desired pulse itself. Equation (A. 18) 
or (A. 29) may be considered as the transmitted signal at 
some reference point (reference plane, antenna terminals, 
etc.). This transmitted signal will propagate into the 
plasma medium and result in a scattered electromagnetic 
field which will then appear at the same reference point. 

One should expect for a frequency dependent medium such as 
a plasma, that the scattered signal will be distorted as 
well as delayed with respect to the transmitted pulse train. 
The effect of the medium will appear in the measured 
reflection coefficient at the reference point for each 
constituent frequency in the signal. This set of reflection 
coefficients can be determined in principle by homodyne 
synchronous (I,Q) detection of the received signal at each 
of the stepped radio frequencies and normalization by the 
magnitude of the transmitted signal. The preceding 
discussion, thus, shows that the response of a time 
independent, but frequency dispersive plasma medium to a 
time domain pulse train signal, can be rigorously emulated by 
a sequence of frequency domain measurements taken one 
frequency at a time. 

If R is the complex reflection coefficient at 
n 

frequency n <*> Q , then the effective received time domain 
signal f^(t) due to equation (A. 18) is 
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f R (t) = 


M+k 

E R 

n=M 


-inw t l (M+k) in t 

n® 0 + 2 E M R n e 

n= -M 


M+k 

= E |R | cos(n<o rt t + 4>J 
n=M 


n 


where 


^ n = tan 


1 - 


Im (R ) 
n 


Re (R 


(A. 30) 


n 


Note that R 


n 


is complex and that we take R =R *. 

-n n 


f R (t) is real. 


Thus , 

As written, equation (A. 30) is the response 
of the medium to a pulse train that resembles the "sine” 
modulated pulse train form of equation (A. 28). A "window" 
sequence W n can be applied at this point to give the 
received windowed signal fp^, which may be written as 


? RW (t) = 2 E 


M+k 


R W e in V + ir; (M+ R ) W Q ~ in< * > 0 t 
n=M n n 2 n=-M n n 


M+k 

= E W R |R n | cos(nw Q t + 4> n ) , (A. 31) 

n=M 


where is the same as in equation (A. 30) and W Note 

n n - n 

again that each R n is measured essentially as a (complex) 

d.c. quantity even though the value of R is that 

n 

appropriate for the response of the medium at frequency no} Q . 
Since we are interested primarily in the time delays 
associated with the received signals in order to measure 
distances to the turning point, we are not really interested 
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in the radio frequency behavior of f^(t) of equation 
(A. 31). Without any loss of desired information, we can set 
M=0 in equation (A. 31) and emulate a baseband pulse train 
which represents the envelope of the received signal. The 
final version of the emulated received signal is a received, 


windowed, baseband signal f 


RWB’ 


written as 


f RWB (t) 


= i E W 0? e in V + \ £ W Re in “o t 

2 n=0 n n 2 n=0 n n 


= E W R cos(n<*> t + $ ) 
** n n 1 n 1 ' o n 

n=0 


(A. 32) 


with 


* = 
n 


tan 


Im(R n ) 
Re (R n ) 


where 


R =R 
n 


M+n 


(A. 33) 


and 


v w 


M+n 


(A. 34) 


The baseband transmitted pulse train may be represented by 
equation (A. 13) as 




k 


= E 

n=0 


W cos nw t 
n o 


(A. 35) 
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Equation (A. 32) will be a distorted and delayed version of 
equation (A. 35). The comparison of the delay between the 
two equations represents the distance measurement. The left 
side of equation (A. 32) can either be obtained by the 
indicated summation for any t or by using the inverse 
discrete Fourier transform on the complex amplitudes IW tfc 
The direct evaluation of equations (A. 32) and (A. 35) can 
produce the values of f° r any value of t, so that a 

smooth curve can be obtained. The discrete Fourier 
transform also tends toward a smooth curve if the number 
of points is increased in the transform by adding zeroes for 
higher frequencies. The zero padded transform does an 
interpolation between points that is just what the Fourier 
series of equation (A. 32) produces directly. 
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TABLE I 


REFLECTION COEFFICIENTS AT FRONT INTERFACE OF 
INHOMOGENEOUS LAYER-LINEAR FROFILE 


74.0000 

74.0317 

74.0635 

74.0952 

74.1270 

74.1587 

74.1905 

74.2222 

74.2540 

74.2857 

74.3175 

74.3492 

74.3810 

74.4127 

74.4444 

74.4762 

74.5079 

74.5397 

74.5714 

74.6032 

74.6349 

74.6667 

74.6984 

74.7302 

74.7619 

74.7937 

74.8254 


Runge 

-Kutta 

Exact 

-Airy 

Real 

Imag 

Real 

Imag 

1.0000 

0.0000 

1.0000 

0.0000 

0.9682 

0.2503 

0.9682 

0.2503 

0.8746 

0.4848 

0.8746 

0.4848 

0.7252 

0.6885 

0.7252 

0.6885 

0.5294 

0.8484 

0.5294 

0.8484 

0.2996 

0.9541 

0.2996 

0.9541 

0.0504 

0.9987 

0.0504 

0.9987 

-0.2023 

0.9793 

-0.2023 

0.9793 

-0.4422 

0.8969 

-0.4422 

0.8969 

-0.6538 

0.7566 

-0.6539 

0.7566 

-0.8235 

0.5673 

-0.8236 

0.5672 

-0.9401 

0.3410 

-0.9401 

0.3408 

-0.9957 

0.9222 

-0.9958 

0.9207 

-0.9867 

-0.1628 

-0.9866 

-0.1630 

-0.9132 

-0.4076 

-0.9131 

-0.4078 

-0.7799 

-0.6259 

-0.7797 

-0.6261 

-0.5953 

-0.8035 

-0.5951 

-0.8036 

-0.3715 

-0.9285 

-0.3712 

-0.9286 

-0.1229 

-0.9924 

-0.1227 

-0.9924 

0.1339 

-0 . 9^10 

0.1341 

-0.9910 

0.3820 

-0.9242 

0.3822 

-0.9241 

0.6050 

-0.7962 

0.6052 

-0.7960 

0.7881 

-0.6155 

0.7883 

-0.6153 

0.9191 

-0.3939 

0.9192 

-0.3937 

0.9893 

-0.1462 

0.9893 

-0.1460 

0.9938 

0.1114 

0.9938 

0.1116 

0.9323 

0.3617 

0.9322 

0.3619 
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TABL E I ( continued) 


junmii 

Runge- 

Real 

Kutta 
I mag 

Exact- 

Real 

Airy 
I mag 

74.8571 

0.8088 

0 . 5081 

0.8087 

0.5882 

74.8889 

0 . 6314 

0.7754 

0.6313 

0.7755 

7 4 . §206 

0.4119 

0.9112 

0.4118 

0.9113 

74.9524 

0.1647 

0.9864 

0.1645 

0.9864 

74.9841 

-0.9377 

0.9956 

-0.9390 

0.9956 i 

75.0159 

-0.3462 

0.9382 

-0.3463 

0.9381 
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Distance ► 


Figure 1 Reflection from a sample plasma 
electron density profile as a funotion of 
distance. The turning point is located at 
h. . 
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(SOT i ITS /die a 


PLASMA ELECTRON OENSIfY PROFILE 



OCHAICE (centimeters) 


Figure 2 Electron density profile predic- 
tion for MRIS . [P. Gnoffo, "Conservation 

Equations and Physical Model for Hypersonic 
Air Flows in Thermal and Chemical Non- 
equilibrium." NASA TP-2867, February 
1989. ] 
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Distance , z 


Fig ure 4 Electron density profile for a 
uniform dielectric (Region II) between free 
apace regions I and III. 
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iTl cz ure 6 Electron density profile for a 
linear layer (Region II) betneen free space 

redone I and III. 
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Figure 7 Sample linear electron density 
profile. 
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REFLECTION COEFFICIENT 



FREQUENCY (GHz) 

74 lo 76 GHz (64 f»equti*:Ut) 


Figure 9 Magnitude o£ the reflection 
coefficient for a 14 cm. thick uniform 


dielectric slab ( £ j 


: _ 0 . 5 ) located 10 cm 


from measurement plane. 
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M2 (tope) 


REFLECTION COEFFICIENT 


X 



Figure 10 Phase of the reflection coef- 
ficient for a 14 cm. thick uniform dielec 
trie elab (« r =-0.5) located 10 cm. from 

measurement plane . 
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MKSIKFLCCEFF. 


TIME DOMAIN RESPONSE 



Figure 11 Time domain response for a 14 
cm. thick uniform dielectric slab 

t e = _o 5) located 10 cm. from measurement 
v r 

plane. 
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TIME DOMAIN (14 CM DIELECTRIC SLAB) 



. Wmtawatf Drt* 

* K«lMr~8**sat ***•• 


TIME OELAY (nanoseconds) 
74 lo 76 GH« (64 rrtqptncUs) 


Figure 12 Propagation delay for a 14 cm. 
thick uniform dielectric slab (( p =-0.5). 

The shift 2tj corresponds to the round-trip 
propagation delay. The delay is measured 
on the negative side of the time axis. 
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Distance 




Figure 13 Reduction of sample electron 
density profile where the critical electron 
density level is N e cb . The turning point 
moves from h^ to h^ and then h-j as the 
profile is reduced. 



Distance 


Figure 14 Modulation of sample electron 
density profile where the critical electron 
density level is N, cr . The turning point 
moves from d^ to d^ as the ripple on the . 
profile changes. 



nWCFtHi REFL CCEPF. 



TIME DOMAIN RESPONSE 



»JS VAfiWICE N PROFI.E 


HUE (nonoseconds) 

74 to 76 GHz (64 frf quencitr) 


Figure 15 Time domain reeponueo for 
reduaed electron density profiles for 74 to 
70 GHz. Shown are responses for 100%, 95% 
and 90% values for the profile. 
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UMIKJLIM 


PLASMA MOVEMENT (. 



Figure 16 Time domain response for a 10% 
reduction in electron density profile. Tli 
95% response 1 b used as a reference. 




TRAHSFCHiRffLCCEFF. 


TIME DOMMN (modulotion 2 m/s ) 



TIME (nanoseconds) 

mod: .03 anvllud* . 1 cm .ovelontfh 74 j 0 ?6 q Mz ( 6 ., | re q uenc | e s) 


Figure 17 Transform reflection coefficient 
for a modulated electron density profile 
for 74 to 76 GHz with modulation amplitude, 
A=.05 and modulation velooity at 2 m/s. 
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TIME 00 MAIN (modulotion 2 m/s) 



TIME (nanoseconds) 

me* .OS ompltud. . 2 cm ••■'•W 74 (o 76 GMz (6< rreque , K j„) 


Figure 18 Time domain response (Power) for 
a modulated electron density profile for 74 
to 76 GHz with modulation amplitude, A=.05 
and modulation velocity at 2 m/a. 
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Figure 10 Transform reflection coefficient for 74 to, 78 ails 
with Modulation amplitude, A=,06 and modulation vAlbolty at 
(a) 2 m/a, ( b) 4 m/a, (o) 8 m/a and (d) 8 m/a. 
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Figure 21 TJme domain reoponoo for HU to H2 ON* with 
modulation amplitude, A=.U5 and modulation voloolty at 
(a) .22 m/n, (b) t.UU m/a, (c) 2.UU m/e and (d) 2. 07 m/e. 
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flcurs 22 Tin# dona In reaponoa 
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<a) 3.33 m/a, W */*• 
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modulation velocity at. 
6.00 m/a and (d) 8.00 m/a. 
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Figure 24 Tine domain responee for 74 to 76 Oils with 
Modulation amplitude, A=.2G and modulation velocity at 
(a) 3.33 m/a, (b).4.0Q m/o, (o) 4.07 m/a and (d) 5.33 m/e. 
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